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equation of state dependence of the radii of the trapped cloud and the collective 
oscillations in the vicinity of a Feshbach resonance. 
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The newly created ultracold trapped Fermi gases with tunable atomic scattering length 
[1-10] in the vicinity of a Feshbach resonance offer the possibility to study highly corre- 
lated many-body systems including the cross-over from the Bardeen-Cooper-Schrieffer (BCS) 
phase to the Bose-Einstein condensate (BEC) of molecules [1,11-15]. 

In this letter we report our investigation of the dynamics of the strongly interacting dilute 
Fermi gas (dilute in the sense that the range of interatomic potential is small compared with 
inter-particle spacing) at zero temperature. As an example of application we consider the 
expansion of the cloud of 6 Li atoms initially confined in an anisotropic harmonic trap, 
study the equation of state dependence of the radii of the trapped cloud and the collective 
oscillations near a broad Feshbach resonance at a magnetic field B = 820 ± 3G [16-18]. 

We consider a Fermi gas comprising a 50-50 mixture of two different states confined in a 
harmonic trap V ex t(r) = (m/2)(ujj_(x 2 + y 2 ) + oo 2 z 2 ). The s-wave scattering length between 
the two fermionic species is negative, a < 0. 

Our starting point is the single equation approach to the time-dependent density-functional 
theory [19]. The basic of this strategy is to construct the following equation 

,9vT> h 2 

ih-z- = V 2 ^ + V ext V + V XC V (1) 

dt 2m v ; 

that yields the same n(r, t) =\ ty(r,t) \ 2 as the original many-fermions system. The dynamics 
of the system is controlled by an effective single-particle potential V xc (f, t). The central 
problem is the approximation for the xc potential. This is in general a nonlocal functional 
of the density that depends on the history of the system (on the behavior of the density at 
times t' < t). 

The simplest approximation is to ignore nonlocality in space and retardation in time. 
This leads to the adiabatic local density approximation 

V xc {f 1 t) = f^} n=nm (2) 

where e is the ground state energy per particle of the homogeneous system and n is the 
density. The right hand side of Eq.(2) is the local density approximation for the ground - 
state xc potential, but it evaluated at the time-dependent density. Notice that in the above 
equation n is the total density of the gas given by the sum of the two spin component. 

The adiabatic local density approximation is a remarkably good approximation if the 
energy gap is much larger than the oscillator energies Hcu z , hu>± [20,21]. It is expected that 



this condition is satisfied for small temperature [20,21]. Here we notice Refs. [22-24] who 
argue that the ground state of the mixture of two species of fermions with different densities 
(mass) contains both a superfluid and a normal Fermi liquid. We do not consider this 
asymmetrical mixture in the letter. 

The ground state energy per particle, e(n), in the low-density regime, k F \ a |<C 1, can 
be calculated using an expansion in power of k F \ a \ 

e(n) = 2E F (^--^-k F | a | +0.055661(A; F | a |) 2 -0.00914(A; F | a |) 3 -0.018604(A; F | a |) 4 +...), 
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(3) 

where E F = an d k F = (S^n) 1 ^ 3 . The expansion (3) is valid for 3-D. For the case 
of dimensions d < 3, it is known that the quantum-mechanical two-body t-matrix vanishes 
at low energies [25]. The first term in Eq.(3) is the Fermi kinetic energy, the second term 
corresponds to the mean- field prediction [26], the next two terms were first considered in 
Refs. [27,28] and Ref.[29], respectively. 

In the a — > — oo limit (the Bertsch many-body problem, quoted in Ref. [30]) e{n) is 
proportional to that of the non-interacting Fermi gas 

e(n) = (l+/3)±*«Z, ( 4 ) 
10 m 

where a universal parameter (3 [9] is negative and | f3 \ < 1 [30-32] . 
We also consider the following approximations for e{n): 

e[n} ~ Ef{ 5 " 1 + (6/(35tt))(11 - 21n2)A; F | a \ h (5) 

and 

/ \ „ ,3 Sikp | a | +5 2 (k F \a\) 2 

e(n) = E F (--2- — - ' — \ ), (6) 

vy v 5 1 + 5 3 k F | a | +5 4 (k F \ a |) 2 w 

where ^ = 0.106103, 5 2 = 0.187515, 5 3 = 2.29188, 5 4 = 1.11616. 

While Eq.(5) [31] reproduces first three terms of expansion (3) in low-density regime 
and approximately valid in unitary limit, (3 = —0.67, Eq.(6) reproduces first four terms of 
expansion (3) in low-density regime and in unitary limit, k F a — > — oo, reproduces exactly 
results of the recent Monte Carlo calculations [32], (3 = —0.56. 

It can be proved [33] that every solution of the equation 

.,dV h* 2T -r t T d(ne(n)) T 
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is a stationary point of an action corresponding to the Lagrangian density 

£o = |(^_*.^ ) + |_|V*r +f (»> + ^„, (8) 
which for \1/ = e l ^ r ' t ^n 1 ^ 2 (f,t) can be rewritten as 

C = h<Pn + ^(Vv 7 ^) 2 + 7^ n W) 2 + < n ) n + V ^tn- (9) 

The only difference from equations holding for bosons [33,34] is given by density depen- 
dence of e(n). We do not consider three-body recombinations, since these processes play an 
important role near p-wave two-body Feshbach resonance [35]. 

Let us first discuss the expansion of the fermionic superfluid in the a — > — oo limit, 
Eq.(4). In the hydrodynamic approximation (neglecting quantum pressure term, ^(Va/ti) 2 , 
in Eq. (9)) the corresponding Euler-Lagrange equation admit the simple scaling solution, 
n{r,t) = n (xi/bi(t)) [20]. We note here that the hydrodynamic behavior of a cold Fermi 
gas [9] is not in general direct experimental evidence for superfluidity [36-38]. 

We take into account the quantum pressure by finding the optimal ground state energy 

rp 3 * , o4/3 3 

[E^v^+V( i +/ 3 ) 1/2Arl/3 n(v /rr ^^) 1/3 ]- ( io ) 



1=1 1=1 

In this case the scaling parameters obey the following equations 

at t = bi(0) = 1 and ^(0) = 0. 

The predictions of Eqs.(ll) for aspect ratio, u z y/l — ^ z b^/ (u±y/l — j±b z (t)), are reported 
in Fig.l and show that the effect of inclusion of the quantum pressure term on the expansion 
of superfluid is about 1%. For the reminder of this letter we will use the hydrodynamic 
approximation. 

Now we consider a general time-dependent harmonic trap, V ext (f, t) = (m/2) Yh=i ^{fyx?, 
and a general e(n). A suitable trial function can be taken as 4>(r, t) = x(0~K m / i^ 1 ) Si=i r )i{p) x 'li 
n(f,t) = no(xi/bi(t))/ Yijbj. With this ansatz, the Hamilton principle, 5 J dt J C d 3 r = 0, 
gives the following equations for the scaling parameters bi 

uj 2 Il n2de (n)/dn] {rl/Ub d 3 r 
4 



where 6j(0) = 1, &j(0) = and Ui = uii(0) fix the initial configuration of the system, 
corresponding to the density n (r). 

The release energy which corresponds to an integral of motion of Eq.(12) is expressed by 

Erei = -J^l / nW)^d\ + / n e(n /Yl bj )d\}, (13) 

and for the case of e(n) oc n 7 

_ 2/i 7j| 1 
^ 5 7 + 2 [ 2^ + n^J J ' 

where /x is the chemical potential. 

Expanding Eqs.(12) around equilibrium {pi = 1) we get in the case of anisotropic trapping 
/u± = A) the following result for the frequency of the radial compression 

mode 

u rad = ^=[4 + 2k + 3A 2 + k\ 2 + ^(4 + 2k + 3A 2 + /tA 2 ) 2 - 4(10 + 6/t)A 2 ] 1/2 , (14) 

where k = / n\d 2 e/ (drifycPr / J nlde/(dn )d 3 r. For an elongated trap, A < 1, we can rewrite 
Eq.(14) as 

uj rad m u±y/4: + 2k. (15) 

Note that Eq.(14) for the case of e(n) oc n 7 was considered in several papers [40]. 

In Fig. 2 we present the calculations of oj ra d using two approximations, Eq.(5) and Eq.(6), 
for the equation of state e(n) (to calculate the ground-state density we have used a highly 
accurate variational approach of Ref.[41]). The curves explicitly show the nonmonotonic 
behavior of uo ra d in the agreement with a schematic interpolation of Ref. [42] . It can be seen 
from Fig. 2 that the difference between two approximations, Eq.(5) and Eq.(6), is less than 
0.7%. 

Our calculated results for the axial cloud size of strongly interacting 6 Li atoms as a 
function of the magnetic field strength B are compared with the recent experimental data 
[10] in Fig. 3. This comparison shows that although both approximations, Eq.(5) and Eq.(6), 
give a reasonable agreement with experimental data, the equation of state from Eq.(6) leads 
to the better description of the experimental curve. We have used the data from Ref. [17] to 
convert a to B. 
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We note here that in general a Feshbach resonance may lead to the density dependence of 
the effective interaction (for bosons cases see, for example, [43,44]). 

In conclusion, we have considered the expansion of the cloud of initially confined 6 Li 
atoms and studied the equation of state dependence of the radii of trapped cloud and col- 
lective oscillations near the broad Feshbach resonance at B — 820 ± 3G. It is shown a non 
monotonic behavior of the radial compression mode frequency and demonstrated that an im- 
portant test of the equation of state can be obtained from the study of the radii of trapped 
cloud in regimes now available experimentally. 

Note added: A recent paper by the Duke University group [45] reports on measurements 
of the radial compression mode frequencies. Our calculations in a very good agreement with 
these experimental data on the BCS side. 

While this work was being prepared for publication, two preprints [46,47] appeared in 
which the authors consider collective modes and the expansion of a trapped superfluid Fermi 
gas in the BCS-BEC crossover. For the negative scattering length case, their results are in 
perfect agreement with ours. 
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Fig. 1. Aspect ratio of the cloud of the N = 7.5 x 10 4 6 Li atoms as a function of time 
after release from the trap (u± = 2ix x 6605Hz, uj z = 2n x 230Hz). The circular dots indicate 
experimental data from the Duke University group [9]. The solid line and the dashed line 
represent theoretical calculations in the unitary limit (a — > — oo) including the quantum 
pressure term and in the hydrodynamic approximation, respectively. 
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Fig. 2. Radial compressional frequency in unit of u>± as a function of the dimensional 
parameter (A^ 1//6 a/a/ lo )~ 1 . In the unitary limit, a — > — oo( 



one expect tu/cjx = ylO/3 ~ 
1.826. The solid line and the dashed line represent the results of theoretical calculations 
using equations of state Eq.(6) and Eq.(5) respectively. 
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Fig. 3. Axial cloud size of strongly interacting 6 Li atoms after normalization to a non- 
interacting Fermi gas with N = 4 x 10 5 atoms as a function of the magnetic field B. The 
trap parameters are u± = 2n x 640Hz, io z = 27i(600B /kG + 32) 1//2 Hz. The solid line, dashed 
line and circular dots represent the results of theoretical calculations using equations of state 
Eq. (6), Eq. (5) and experimental data from the Innsbruck group [10], respectively. 
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